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Abstract. Topics related to the construction, phenomenological determination, and effects 
of the effective three-body forces within the traditional nuclear shell model approach are 
discussed. The manifestations of the three-body forces in realistic nuclei in the O/7/2 and 
IsOd shell model valence spaces are explored. 



1 Introduction 

In this work we investigate the role that three-body forces play within the nuclear 
shell model (SM) approach. Establishment of the effective interaction parameters, 
study of hierarchy in strength from single-particle (s.p.) to two-body, three-body, 
and beyond, manifestations in energy spectra and transitions rates, comparison with 
different traditional SM calculations, and overall assessment for the need of beyond- 
two-body SM are the topics for this discussion. Previous works in this direction 
have shown an improved description of nuclear spectra [1-6] and the significance 
of three-body monopole renormalizations [7]. 

The effective interaction Hamiltonian of rank fc is a sum 



H, = ^H(-\ where il(")=^^Ff ^ (a) rl"2(/3), 



L 

^^^'^^TT(n) TT(n) T An) / _ n\ ' 

n=l Q/3 L M=-L 

(1) 

is the n-body rotationally invariant component of the interaction. The 71 -particle cre- 
ation operators t'^^^ [a) with the total angular momentum L and magnetic projec- 
tion M are normalized (0|r^?]^, (a') {^)\^) = ^aa'SLL'^MM' and expressed 

through the s.p. creation operators as T^Jj (a) — J2i2 n ^i2^.n((^) o^iO't- ■ - ^lii 
where index 1 labels a s.p. state. The choice of coefficients CiJ^ j^{a) that de- 
fines a full set of orthogonal operators T^'^(a) is generally not unique. For nu- 
merical work it is most convenient to use a full set of orthogonal eigenstates 

|n; LMa) = T^^j (a)|0) of some n-particle system [5]. In the m-scheme SM we 
generate states only for a particular value of the total magnetic projection A/, the 
remaining states are obtained by the raising and lowering angular momentum opera- 
tors. It is possible [5], to select a single reference two-body Hamiltonian which then 

can be used to define all many-body operators T^"^ (a) for n > 2. The traditional 
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SM Hamiltonian is H2 = H^^^ + H^"^") , where the two-body operators T^lj (a) are 
determined with the help of the Clebsch-Gordan coefficients. 

2 Manifestation of three-body forces in f 7/2 -shell nuclei 

As a first example we present here a study of a single-j O/7 /2 shell, related discus- 
sion may be found in Ref. [6]. We consider two types of systems = 28 isotones 
starting from ^^Ca with protons filling the O/7 /2 shell and the Z=20, ^°^^^Ca iso- 
topes with valence neutrons. The states in these systems that are identified by exper- 
iments with the /■j/2 valence space are listed in Tab. [T] The shell is unique 
because of symmetries associated with the quasispin and particle-hole conjuga- 
tion [1,8-11]. These symmetries are violated if interaction is beyond the two-body. 

2.1 Particle-hole symmetry 

The violation of the particle-hole symmetry is due to monopole terms that are 
non-Uner in the particle-number density. These terms in the Hamiltonian appear 
from three-body and higher rank interactions [7]. For a single-j and a standard 
two-body SM the symmetry is exact and it makes the spectra of N and N — 
fl — N particle systems identical, apart from a constant shift in energy, here 
— 2] + 1. The particle-hole conjugation operator C that acts on a s.p. state 
as a]^ = Caj„C~^ = (— l)-'^™aj_„,, transforms an arbitrary ri-body interac- 
tion into itself plus some Hamiltonian of a lower interaction-rank namely 
J?'^"^ = (— + -ff.n-i- The n = 1 case represents a particles to holes trans- 
formation N = —N + Q. For the n = 2 it leads to a monopole shift 

fj(2) ^ jj(2) ^ - 2N)M, M = ^(2L + l)V^\ (2) 

Within a single-j one-body Hamiltonian is a constant of motion, being always pro- 
portional to N. Thus, following Eq. the two-body interaction is identical for par- 
ticles and holes, apart from some constant-of-motion term. The interaction of rank 

3 and higher violate this symmetry making excitation spectra of N and N = fi — N 
particle systems different. The experimental data in Tab. [T] shows the particle-hole 
symmetry violations, for example the excitation energies of = 2 states in = 2 
system are systematically higher then those in the 6-particle case, indicating a re- 
duced ground state binding. Using this information a monopole component of the 
three-body force can be extracted from the differences in excitation energies be- 
tween particle and hole systems, see Fig.[T|and discussion below. 

2.2 Seniority 

The j — 7/2 is the largest single-j shell for which the number of unpaired nucleons 
V, the seniority, is an integral of motion for any one- and two-body interaction [8, 
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Table 1. States in f-j/2 valence space with spin and seniority listed in the first and second 
columns. The * denotes seniority mixed states in 38/7/2- Following are columns with data 
for N = 28 isotones and Z — 20 isotopes. Three columns for each type of valence particles 
list name and excitation energy, experimental binding energy, and energy from the three-body 
SM calculation discussed in the text. All data is in units of MeV. 
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(2) (2)t 

12]. Formally, the pair operators Tqo ' -^oo ' '^^^ particle number N form an SU(2) 
rotational group, which because of its analogy to angular momentum is referred to 
as quasispin. The relation is established by the operators 

_iv n n V 
^^Y^T' T"2 

with £(£ + 1) being an eigenvalue of the quasispin vector squared and Cz its mag- 
netic projection. 

For a spectrum, the invariance under seniority sets relations between states of the 
same L but different projection C^- For example, the excitation energies of = 2 
states from the = 0, 0+ ground state are identical in all even-particle systems. 
Using Wigner-Eckart theorem a full set of relations can be established, see for ex- 
ample sec IIIB in ref. [11] or Ref. [13]. The invariance under quasispin rotations 
allows to classify operators in close analogy to the usual rotations. The s.p. opera- 
tors associated with the particle transfer reactions carry £ = 1/2 and thus permit 
seniority change Av = 1. The reactions ^^Wi^He, d)^^Cr and ''■^Ca(d,p)^'*Ca in- 
dicate seniority mixing as ly — A final states are populated [14, 15]. The one-body 
multipole operators are quasispin scalars for odd angular momentum, and quasispin 
vectors for even. Thus, the Ml electromagnetic transitions are given by the qua- 
siscalar operators that do not change quasispin. The E2 operator is a quasivector 
In the mid-shell for ^'^Cr and '^^Ca, where N = f2/2 and ^ the E2 transi- 
tions between states of the same seniority are forbidden. Seniority can be used to 

(n) 

classify the many-body operators T^^'j and interaction parameters. The three-body 
interactions mix seniorities, one exception is the interaction between 1^ = 1 nucleon 

(3) 

triplets given by the strength V^,^ . 



2.3 Parameter fit and evidence for three-body forces 

To obtain the parameters of the effective Hamiltonian with the three-body forces we 
conduct a full least-square p ~ll parameter fit to d =31 data-point (27 in the case 
of Ca isotopes). The procedure is similar to a two-body fit outlined for this model 
space in sec. 3.2 of Ref. [16]. Schematically E = AV where E is a set of 31 ener- 
gies, V is a list of 1 1 interaction parameters and A is 3 1 by 11 matrix created from 
the linear form of the Hamiltonian operator Eq. ([TJ. Due to the seniority mixing A 
depends on the eigenstates, which in turn are determined by the interactions V; thus 
the overall fitting procedure is iterative [17]. In this example seniority mixing occurs 
only in 4 states of the = 4 system, in consequence most of the matrix elements 
of A are constants. Using the set of experimental data in Tab. [T] denoted here as 
Eex, we determine V = (A-^A)^^ A-^Ee^^, where T and —1 superscripts indicate 
transposed and inverted matrices. The obtained interaction parameters are used to 
update non-constant elements of A. The procedure is iterated several times so that 
the interaction-dependent components of A converge. In Tab.|2]the resulting param- 
eters are listed for the = 28 proton system, and for the neutron Z = 20 system. 
The two columns in each case correspond to fits without (left) and with (right) the 
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three-body forces. The root-mean-square deviation (RMS) \Eex — AV| / ^/d is given 
for each fit. The confidence Umits given in brackets are inferred from the variances 
for each fit parameter 
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Table 2. Interaction parameters of 2B/7/2 and 3B/7/2 SM Hamiltonians determined with the 
least-square fit are given in keV. 



The reduction of the RMS deviation, for example for Z — 28 isotones it drops 
from 120 keV to about 80 keV, is not the only evidence in support of the three- 
body forces. The fit parameters are stable within quoted error-bars even if some 
questionable data-points are removed. The energies from the three-body fit listed 
in Tab. [T] are comparable or even better than the results from many two-body SM 
calculations in the expanded model space [18, 19]. However, such comparisons are 
difficult since different models have different number of parameters and were fit to 
different sets of nuclei. 

In Tab. [3] we discuss the renormalization of pairing by considering a minimal 
fit limited to the ground states and a single three-body term. The fit is similar to 
Ref. [8], but has a seniority conserving three-body force given by the v ^ \ triplet 
operator T^J^ ~ ajm^oo "^^^^ the strength Vy!^. This interaction is equivalent to 
a density-dependent pairing force [20]. In a single-j shell the renormalization of 
pairing by a particle-number dependent strength 

vr = vr + ^^^f (3) 
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allows for an exact treatment of the three-body term. The ground state energies with 
= or 1 are 

E^eN+^^j-^(^{[2-N-iy)^^ + {N-2 + i^)M'^, (4) 

which is a usual expression [8, 11], but includes a renormalized pairing strength 
denoted with prime. The results from the minimal fit are shown in Tab. |3] they are 
consistent with the full fit in Tab. |2] 
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Table 3. Interaction parameters for the minimal f-j/2 SM determined with the linear least- 
squared fit of 8 binding energies. In brackets the variances for each parameter are shown. The 
two columns for isotopes and isotones are fits without and with the three-body term. 



In Fig. [T]we give a cumulative picture showing the Vy}^ term found with differ- 
ent methods. As discussed above, due to the particle-hole symmetry and seniority 
conservation, excitation energies of = 2: 2+, 4+, and 6+ states in = 2, 4, and 
6 -particle systems should be identical. The Vy!^ coefficient can be found assuming 
that it is responsible for most of the mass difference. For example, the difference 
in excitation energies of these states between ^"Ti and ^^Cr equals to diVyl^/i. The 

independent result on Vy^ inferred from these observations, the binding energy fit, 
and the fit to all states in the = 28 isotones with 4 parameters are summarized in 
Fig.[T] The point that corresponds to the 4+ state in ^^Cr in Fig.[T|is not in agreement 
with the rest of the data, it demonstrates the seniority mixing discussed below. 

It follows from Tabs. |2] and [3] and Fig.[T]that within the error-bars the three-body 
interaction is isospin invariant; it is the same for proton and neutron valence spaces. 



2.4 Seniority mixing in ®^Cr 

The mid-shell case of ^^Cr, see Fig|2] is interesting to discuss. Here, in addition to 
2B/7/2 and 3B/7/2 interactions from Tab.|2]we perform a large scale SM calcula- 
tion 2B/7/2P (includes pi/2 and P3/2) and 2B/p (entire /p-shell, truncated to 10'' 
projected m-scheme states) using FPBP two-body SM Hamiltonian [21]. Similar 
results in a more restrictive valence space can be found in Ref. [18]. 

The level repulsion between neighboring 4j'" and 4^ states is generated by the 
seniority mixing, the observed energy difference of 400 keV is not reproduced by the 
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Figure 1. Cumulative data on Vy^ seniority v = 1 effective three-body force in A'' = 48 
isotones, left, and Z = 40 isotopes to the right. The point labeled as BE comes from a fit to 
8 binding energies in Tab[3]and includes a fitting error-bars. The point labeled as "Full Fit" 
corresponds to a fit of all 3 1 levels in TV = 28 isotones with 6 parameters for s.p. energy, two- 



body force and Vyl . Other individual points correspond to extraction of Vy^ from excitation 
energy, always compared to A'^ = 2 system (^"Ti or '''^Ca) 
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Figure 2. Spectrum of 52Cr. 
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2B/7/2 (84 keV) model. As seen in Fig.|2]the discrepancy remains in the extended 
two-body model 2B/7 /2P (200 keV). Although, the full 2B fp model reproduces the 
splitting, the excessive intruder admixtures over-bind the ground state and effec- 
tively push all states up in excitation energy. The 3B/7/2 model is in good agree- 
ment with experiment; its predictions for the seniority mixing are I'iii ) = 2.82 
and 1^(4^) = 2.71, as calculated from the expectation value of the pair operator 

(^oo^^Too^) = (N -:^){2j + 3- N - i/)/(4j + 2). The 2+ state is relatively pure 
1/(2+) = 2.006. 

The seniority mixing violates quasispin selection rules [10, 15,22-24] which in 
the past have been explained by the two-body models beyond the single-j [9, 10, 
14, 18, 25], however such models not always describe all of the features observed in 
experiment. In particular, to explain electromagnetic transitions sizable variations 
of effective charges are needed [26] and the particle transfer spectroscopic factors 
do not show large amount of strength outside the /7/2 valence space [22]. In Tab. 
^B{E2) transitions rates from all models are compared to experiment. To make a 
fair comparison the combination of the nuclear radial overlap and effective charge 
is normalized using observed E2 rate for the transition 2i — > Oi in the 2B/7/2, 
2B/7/2P, and 2B fp models. The parameter for the 3B/7/2 model is identical to the 
one used in the 2B/7/ 2- The small difference in 2 1 Oi i3(£'2) between the 3B/7/2 
and 2B/7/2 models is a result of the 1/ — 4: admixture in the 2^ state. The strong 
V = 4: and ly = 2 seniority mixing between 4+ and 4^ states impacts forbidden 
transitions; for example, E2 transitions 42 2i and 61 A2 become allowed. 
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57.8 
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27.6 


51.1 


30±10<^* 



Table 4. B(E2) transition summary on ^^Cr expressed in units e^fm'*. The data is taken from 
[27]. '''in the 2B /7/2P and 2B/p models we use 0.5(neutron) and 1.5(proton) effective 
charges, the overall radial scaling is fixed by the B(E2,2i — > Oi). ^^'The life-time error-bars 
were used. '^^^ There are conflicting results on life-time; we use DSAM (HI, X717) data from 
Ref. [27], which is consistent with [26]. 

The proton removal spectroscopic factors in Tab.|5]show a similar picture, where 
the seniority mixing has a strong impact on transitions. In support of the three-body 
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forces as a source of the mixing it was argued in Ref. [22] that the sum of spectro- 
scopic factors for 4+ states is close to 4/3 which is consistent with the observation 
in Ref. [22] and does not support the expanded valence space where spectroscopic 
factors are reduced due to fragmentation of the single-particle strength. 
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Table 5. Proton removal spectroscopic factors. The experimental data is taken from 
^^V(^He,d)^^Cr reaction [22]. Within error-bars this data is consistent with results [27]. 



3 Three-body forces in oxygen isotopes 

The above single-j example is remarkable due to its transparency and simplicity. 
The general SM case, however, is complicated by an enormously large number of 
parameters and thus difficulty of the fit [2^,7] . Selecting dynamically relevant com- 
ponents of the many-body forces requires an in-depth microscopic understanding 
of their origin. Establishment of the physically relevant set of the operator basis 

^LM ('^) important start. As discussed in the introduction, for n > 2 the index 
a must include an additional information about the coupling scheme, the choice of 
which is not unique. Previous ideas on selecting the best set of triplet operators in- 
clude a possibility of using the v — 1 operators for each single-particle level [20]. 
For i = 7/2 the three-body force associated with this operator, discussed in Fig. 
[T] is indeed a dominating component in binding. However, it is not clear if such 
construction, built upon s.p. levels, is the best choice in a general case given renor- 
malization of the s.p. effective degrees of freedom by the two-body interaction. The 
two-body Hamiltonian of the pairing type would, for instance, suggest the use of 
quasiparticles. 

Perusing this idea we propose an alternative approach which assumes a hierar- 
chy of forces, where higher rank components of the Hamiltonian are perturbative, 
and the operator basis are selected using the many-body dynamics. 

Consider Hn-i, n > 3 Hamiltonian to be determined by some procedure. While 
building a higher rank forces !!„ = H„-i + H^"\ we assume if^") to be per- 
turbatively small. Thus, within the lowest order perturbation theory the n-particle 
wave-functions of Hn-i and iJ„ are the same and can be found by diagonalizing 

Hn-i\n]LMa) = £'„;L(a) \n]LMa). 
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We use these eigenstates to define a full set of n-body operators Tj^"^j{a) as 
|n; LMa) = which we view as the most relevant basis. With a per- 

turbative nature in mind the term i? is diagonal in these basis, V^"' (a/?) = 
if a 7^ /?. Thus, the number of parameters is reduced. Further steps can be taken 
to discuss the significance of the diagonal parameters. When pairing is impor- 
tant one can take only those states (basis operators) that correspond to the lowest 
quasiparticle excitations. Experimental data can be used for guidance. For exam- 
ple, if the n-particle states are known and identified experimentally to have energies 
j^iexp) j.^. ^ direct fit can be done by setting the corresponding n-body interac- 
tion parameters as V^"^ (a, a) = E^^.^^ (a) — En-L (a) , so that the new Hamiltonian 
reproduces exactly the experimental energies. 

There are some issues to stress. Certainly, the transition from n = 1 to n = 2, 
is not a subject to this approach. One has to have a starting SM Hamiltonian H^g}^ 
determined from G-matrix techniques or by other methods, see [21] and references 
therein. It is possible to rewrite the two-body Hamiltonian as a diagonal structure 
by introducing new pair operators, this is useful for perturbative adjustments of the 
two-body interactions. 

The two-body SM Hamiltonian can be used as a primary component of inter- 
action, defining many-body operators, and treating all higher rank forces as per- 
turbations. It is important for this approach to stay within the perturbation theory. 
Departing a perturbative form, it is feasible with this construction to create a Hamil- 
tonian that exactly reproduces energies of all states within a given valence space, 
tests show that in this case the many-body forces have an inverse hierarchy with 
higher rank ones giving a bigger contribution. 

It is an established practice in the SM approach to include a mass dependence 
of the two-body forces. For a short range delta-type interaction the radial overlap 
integrals scale as R~^^^, where R is the radius of the nucleus. Thus, in terms of 
the mass-number A the two-body interaction iJ^^^ ~ A~^/'^. At the opposite ex- 
treme the long-range Coulomb leads to an A~^^^ scaling. The fits to experimental 
data lead to a compromising middle value A~'^^^ [17]. The many-body forces are 
expected to be short range, requiring all participating particles to be localized. The 
resulting scaling that follows from this argument is 

i^(n)^^(l-n)/2^S0 JjW (A) = f ^i£±^ j H^''\A,), (5) 

where Ac is the mass of the core. At this stage it is not clear if this argument is valid 
and if scaling should be included. 

As a demonstration we discuss here a 3-body force in the case of oxygen iso- 
topes. For the two-body interaction Hamiltonian we take a USD shell model [28]. 
The total number of triplet operators, not counting magnetic projections, is 37 which 
in a general 3-body interaction Hamiltonian gives a large number of parameters. Ex- 
amination of experimental data for ^^O and results from the different shell model 
Hamiltonians USD, USDA and USDB [17] show a rather systematic difference; in 
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particular for the lowest 5/2 j^, 3/2 and l/2'l states. These are one quasiparti- 

cle excitations. Thereby, we define the corresponding triplet operators with 
j = 5/2, 3/2, and 1/2 from the three-particle eigenstates of the USD Hamiltonian; 
and specify the three-body interaction in the diagonal form 

i/(3) = V V^-^^^ V T^"'' t'") (6) 

j=5/2,3/2,l/2 m=-j 

For Fig.[3]we fit the three parameters in Eq. ^ to the ground states in even systems 
and to the three lowest states with one unpaired particle in the odd systems for mass 
A — 19 io 22 oxygen isotopes. The values from the best fit are V^'j\ = 45 keV, 

^3/2 = -179 keV, and V^f^ = -231 keV. 

The improvement in the spectrum, seen in Fig. [3] is significant. Certainly, this 
first study is to be continued, there is a possibility to examine more interaction terms, 
discuss scaling of the matrix elements, and to consider fitting all parameters for one- 
, two-, and three-body components together. Modifying perturbatively the two-body 
part should not invalidate the quality of the USD-defined three-body basis. 
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Figure 3. Spectrum of ^'^~'^'^Q isotopes. For every nucleus, experimentally observed states 
are compared with the spectrum that includes three-body forces and with the two-body USD 
SM Hamiltonian, from left to right, as identified at the bottom. 



4 Conclusion 



Dealing with many-body forces, understanding their origins, structure, and hier- 
archy of renormalizations is an important component for a successful solution of 
a many-body problem. This presentations aims to continue the discussion in Ref. 
[1-7] related to the phenomenological three-body forces within the context of the 
nuclear shell model approach. The study of nuclei in the O/7 /2 shell shows evidence 
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of such forces through an overall fit to data with full examination of uncertainties, 
via examination of binding energies and associated differences in excitation spectra, 
and with an in-depth analysis of violations of symmetries in the structure of wave 
functions. 

The general SM problem with many -body forces is complicated by a large num- 
ber of parameters, the absence of a good microscopic approach, difficulties in fits 
and questions related to renormalizations of strengths. These issues are discussed 
and some methods for dealing with them are proposed. In particular, in analogy to 
a Hartree-Fock procedure where single-particle states are defined in the way to best 
represent the dynamics of the system, we propose here methods to identify the most 
relevant many-body operators. These techniques are demonstrated using a chain of 
oxygen isotopes. 

Support from the U. S. Department of Energy, grant DE-FG02-92ER40750 is 
acknowledged. 
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